DYNAMICS AND ABSTRACT COMPUTABILITY: COMPUTING 

INVARIANT MEASURES. 



STEFANO GALATOLO, MATHIEU HOYRUP AND CRISTOBAL ROJAS 

Abstract. We consider the question of computing invariant measures from an abstract 
point of view. We work in a general framework (computable metric spaces, computable 
measures and functions) where this problem can be posed precisely. We consider invariant 
measures as fixed points of the transfer operator and give general conditions under which 
the transfer operator is (sufficiently) computable. In this case, a general result ensures 
the computability of isolated fixed points and hence invariant measures (in given classes of 
"regular" measures). This implies the computability of many SRB measures. 

On the other hand, not all computable dynamical systems have a computable invariant 
measure. We exhibit two interesting examples of computable dynamics, one having an SRB 
measure which is not computable and another having no computable invariant measure at 
all, showing some subtlety in this kind of problems. 



Contents 

1. Introduction 2 
1.1. Plan of the paper 4 

2. Preliminaries on algorithmic theory 5 

2.1. Analysis and computation 5 

2.2. Background from recursion theory 5 

2.3. From N to countable sets 6 

2.4. Computability of reals 6 

2.5. Computable metric spaces 7 

2.6. Recursively compact sets: approximation from above l9 

2.7. Recursively precompact 11 

2.8. Recursively closed sets: approximable from below 12 

3. Computable measures 13 

4. Dynamical systems, statistical behavior, invariant measures 15 

4.1. The transfer operator 16 

4.2. Computing invariant "regular" measures 16 

4.3. Unbounded densities 18 

5. Computable systems having not computable invariant measures 18 
5.1. A computable system having no computable invariant measure 19 
References 21 



Date: March 12, 2009. 

1 



1. Introduction 



An important fact motivating the study of the statistical properties of dynamical systems 
is that the pointwise long time prediction of a chaotic system is not possible, while the 
estimation or forecasting of averages and other long time statistical properties is sometimes 
possible. This often corresponds in mathematical terms to computing invariant measures, 
or estimating some of their properties. 

Giving a precise meaning to the computation of a continuous object like a measure is not 
a completely obvious task and involves the definition of effective versions of several concepts 
from mathematical analysis. 

Our approach will be mainly based on the concept of computable metric space. To give 
a first example, let us consider the set M of real numbers. Beyond Q there are many other 
real numbers that can be handled by algorithms: vr or \pl for instance can be approximated 
at any given 'precision (with rational numbers) by an algorithm. Hence these numbers can 
be identified with the algorithm which is able to calculate them (more precisely, with the 
string representing the program which approximates it). This set of points is called the set 
of computable real numbers and was introduced in the famous paper |T36j . 

This kind of construction can then be generalized to many other metric spaces, considering 
a dense countable set that plays the same role as the rationals in the above example. Then, 
computable or recursive counterparts of many mathematical notions can be defined, and 
rigorous statements about the algorithmic approximation of abstract objects can be made, 
also obtaining algorithmic versions of many classical theorems (see section 2). In particular, 
this general approach also gives the possibility to treat in a simple way measures spaces, 
define computable measures and computable functions between measure spaces (transfer 
operators), which will be the main theme of this paper. 

The paper is devoted to the problem of computation of invariant measures in discrete 
time dynamical systems. By discrete time dynamical system we mean a system {X, T) 
were X is a metric space and T : X —>■ X is a. Borel measurable transformation. Here 
an invariant measure is a Borel measure fi on X such that for each measurable set A it 
holds fi{A) = fi{T~^{A)). Such measures contain information on the statistical behavior 
of the system {X, T) and on the possible behavior of averages of observables along typical 
trajectories of the system. The map T moreover induces a function Lj- : PM{X) — > PM{X), 
where PM{X) is the set of Borel probability measures over X and will be endowed with a 
suitable metric (for details see section [3]). Lt is called the transfer operator associated to T 
(basic results about this are reminded in section H]). 

Before entering into details about the computation of measures and invariant measures 
in particular, we remark that whatever we mean by "approximating a measure by an al- 
gorithm", there are only countably many "measure approximating algorithms" whereas, in 
general, a dynamical system may have uncountably many invariant measures (usually an 
infinite dimensional set). So, most of them will not be algorithmically describable. This is 
not a serious problem because we can put our attention on the most "meaningful" ones. An 
important part of the theory of dynamical systems is indeed devoted to the understanding 
of ''physically" relevant invariant measures, among these, SRB measures play an important 
rolcl. These measures are good candidates to be computed. The existence and uniqueness 

"'^Informally speaking, these are measures which represent the asymptotic statistical behavior of "many" 
(positive Lebesgue measure) initial conditions, see section |4] 
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of SRB measures is a widely studied problem (see |Y02] ) . which has been solved for some 
important classes of dynamical systems. 

Let us precise the concept of computable measure. As mentioned before, the framework 
of computable analysis can be applied to abstract spaces as the space PM{X). A measure 
/i is then computable if it is a computable point of that measure space. In this case there 
is an algorithm such that, for each rational e given as input, outputs a "finite" measure (a 
finite rational convex combination of Dirac measures supported on "rational" points) which 
is £-close to fi. 

In the literature, there are several works dealing with the problem of approximating in- 
variant measures, more or less informally from the algorithmic point of view (see e.g. [LOlj . 
|H95j . |KMY98] . |PJ99] . |Din93l IDin94j ). In these works the main technique consists in an 
adequate discretization of the problem. More precisely, in several of the above works the 
transfer operator associated to the dynamics (see subsection 14. ip is approximated by a finite 
dimensional one and the problem is reduced to the computation of the corresponding relevant 
eigenvectors (some convergence result then validates the quality of the approximation). 

Another strategy to face the problem of computation of invariant measures consist in 
following the way the measure /x can be constructed and check that each step can be realized 
in an effective way. In some interesting examples we can obtain the SRB measure as limit 
of iterates of the Lesbegue measure /i = lim„^oo L'^ijn) where m is the Lesbegue measure 
and Lt is the transfer operator associated to T. To prove computability of /i the main point 
is to recursively estimate the speed of convergence to the limit. This sometimes can be 
done using the decay of correlations (see |GHR09b] where computability of SRB measures in 
uniformly hyperbolic systems is proved in this way, see |GP09j for general relations between 
convergence of measures and decay of correlations with a point of view similar to the one of 
the present paper). 

Let us illustrate the main results of the paper. Informally speaking, a function T : X X 
is said to be computable if its behavior can be described by some algorithm (for the precise 
definitions see sections 12.51 and I3.0.2p . In this case the pair (X, T) is called a computable 
dynamical system. In this context, the general problem we are facing can be stated in the 
following terms: 

Problem 1. a) Given a computable dynamical system {X,T) does the set of invariant 
measures contain computable points ? 
b) Can they be found in an algorithmic way, starting from the description of the system? 

We will see that, in general, even the above question a) does not always have a positive an- 
swer. However, in many interesting situations, both of the above problems can be positively 
solved. 

We will take a general point of view finding the interesting invariant measure as a fixed 
point of the transfer operator, giving general conditions ensuring its computability. The 
following theorem will be the main tool (see Thm. I4.2.ip . 

Theorem A Let X be a computable metric space and T a function which is computable on 
X\D. Let us consider the dynamical system {X,T). Suppose there is a recursively compact 
set of probability measures V C PM{X) such that for every fi & V , ^{D) = holds. Then 
every invariant measure isolated (in the weak topology) in V is computable. 
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The precise meaning of computability on X \ D will be given in section 12.51 however the 
intuitive meaning of the above proposition is that: if the function T is computable outside 
some singular set D (the discontinuity set for example) and we look for invariant measures 
in a set V of measures giving no weight to the set D (some class of regular measures e.g.) 
and in the set V there is a unique invariant measure, then this measure can be computed. 

This will give as a consequence that the SRB measure is computable in many examples of 
computable systems (uniquely ergodic systems, piecewise expanding maps in one dimensions, 
systems having an indifferent fixed point and many other systems having an unique absolutely 
continuous invariant measure, see Theorem 13.0.21 and Prop. 14.2.21 ). 

Observe that any object which is "computable" in some way (as T, V, fi in the theorem) 
admits a finite description (a finite program). Theorem A is actually uniform: there is 
a single algorithm which takes finite descriptions of T and V and which, as soon as the 
hypothesis in Theorem A are satisfied and is a unique invariant measure in V, outputs a 
finite description of fi (see remark |4 . 2 . 1 1 and the above item b) of Problem [T]). Observe that 
the algorithm cannot decide whether the hypotheses are satisfied or not, but computes the 
measure whenever they are fulfilled. 

After such general statements, one could conjecture that, in computable dynamical sys- 
tems, SRB measures are always computable. This is not true, and reveals some subtlety 
about the general problem of computing an invariant measure. In section [5] we will see that: 

Examples There exists a computable dynamical system having no computable measure at 
all. Moreover, there exists a computable dynamical system on the unit interval having an 
SRB measure which is not computable. 

The interest of the second example comes from the fact that any computable map of the 
interval must have some computable invariant measure. The example shows that important 
invariant measures can still be missed. 

To further motivate these results, we finally remark that from a technical point of view, 
computability of the considered measure is a requirement in several results about relations 
between computation, probability, randomness and pseudo-randomness (see e.g. |LM08j . 
|GHR09a] . |GHR09bj . pHR09cj l 

1.1. Plan of the paper. In section 2 we give a compact and self contained introduction 
to the prerequisites about computable analysis which are necessary to work with dynamical 
systems on metric spaces. In this section we also prove some general statements about 
solutions of equations on metric spaces which will be used to "find" the interesting invariant 
measures as fixed points of the transfer operator (Theorem, I2.6.3p . 

In section 3 we develop the computable treatment of the space of probability measures on 
a given (computable) metric space. Some results of these initial sections are new and should 
be of independent interest. Their usefulness is demonstrated in the next sections. 

In section 4 we start considering dynamical systems. A direct application of the results 
of the previous sections allow us to establish general assumptions under which the transfer 
operator is computable (on a suitable subset. Theorem 14.1.11) . 

We then use the framework and tools introduced before to face Problem [TJ We prove 
Theorem A above (which also becomes a simple application of previous results) and show 
how to apply it in order to prove the computability of many interesting invariant measures. 
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In section 6 we construct the two counter-examples already announced. 

2. Preliminaries on algorithmic theory 

2.1. Analysis and computation. A way to approach several problems from mathematical 
analysis by computational tools is to approximate the "infinite" mathematical objects (ele- 
ments of non countable sets, as real numbers or a functions ) involved in the problem by some 
algorithm which constructs an approximating sequence of "finite" objects (rational numbers, 
polynomials with rational coefficients) which are "treatable" by the computer. Usually, the 
algorithm has to manipulate and decide questions about the various mathematical objects 
involved, and convergence results should be provided in order to choose the suitable level 
of accuracy for the finite approximation. The actual implementation of the algorithm and 
the various decisions are, in most cases, subjected to round off errors which can produce 
additional approximation errors, wrong decisions or undecidable situations if the error is 
not considered rigorously (how to decide x > y when x = yl). Sometimes, estimates (for 
these errors) can be obtained under suitable conditions, but this is in general a further and 
often nontrivial task (see e.g. |Bla94j ). In this paper we will work in a framework where 
the algorithmic abilities of the computer to represent and manipulate infinite mathematical 
objects are taken into account from the beginning. In this framework (often referred to as 
Computable Analysis) one can rigorously determine which objects can be algorithmically 
approximated at any given accuracy (these will be called computable objects), and which 
can not. 

Here, the word computable is used, but may be adapted to each particular situation: for 
instance, "computable" functions from N to N are called recursive functions, "computable" 
subsets of N are called r.e sets, etc. 

2.2. Background from recursion theory. The starting point of recursion theory was to 
give a mathematical definition making precise the intuitive notions of algorithmic or effective 
procedure on symbolic objects. Every mathematician has a more or less clear intuition of 
what can be computed by algorithms: the multiplication of natural numbers, the formal 
derivation of polynomials are simple examples. 

Several very different formalizations have been independently proposed (by Post, Church, 
Kleene, Turing, Markov. . . ) in the 30's, and have proved to be equivalent: they compute 
the same functions from N to N. This class of functions is now called the class of recursive 
functions. As an algorithm is allowed to run forever on an input, these functions may be 
partial, i.e. not defined everywhere. The domain of a recursive function is the set of inputs 
on which the algorithm eventually halts. A recursive function whose domain is N is said to 
be total. For formal definitions see for example |Rog87 . 

With this intuitive description it is more or less clear that there exists an effective procedure 
to enumerate the class of all partial recursive functions, associating to each of them its Godel 
number. Hence there exists a universal recursive function : N — N satisfying for all 
e, n G N, (p^{{e,n)) = (^^{n) where e is the Godel number of v?e ^^id (■, ■) : — > N is some 
recursive bijection. 

The notion of recursive function induces directly an important computability notion on 
the class of subsets of N: a set of natural numbers is said to be recursively enumerable 
{r.e for short) if it is the range of some partial recursive function. That is if there exists 
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an algorithm listing (or enumerating) the set. We denote by Ee the r.e set associated to ip^, 
namely: Ee = range(v9g) = {(Pu{{e, n)) : n G N}, where ip^ is the universal recursive function. 

Let {Ei)i^^ be a family of r.e subsets of N. We say that Ei is r.e uniformly in i if there 
is a single recursive function ip such that E^ = {ip{{i,n)) : n G N}. Taking ip = ip^ the 
universal recursive function yields an enumeration {Ei)i^^ of all the r.e subsets of N, such 
that Ei is r.e uniformly in i. 

More generally, once a computability notion has been defined for some class of objects in 
the following form: 

An object x is computable if there is a (partial or total) recursive function ip which 

computes x in some sense. 

A uniform version will be implicitly defined and intensively used: 

Objects from a family (xj)jgN of X are uniformly computable if there is a single (total or 
partial) recursive function ip such that (/^((i, .)) : N — > N computes for each i. 



2.3. Prom M to countable sets. Strictly speaking, recursive functions only work on natural 
numbers, but this can be extended to the objects (thought of as "finite" objects) of any 
countable set, once a numbering of its elements has been chosen. 

Definition 2.3.1. A numbered set O is a countable set together with a surjection uq '■ 
N ^ O called the numbering. We write o„ for //(n) and call n a name of o„. j 

The set Q of rational numbers can be injectively numbered Q = {qo, gi, . . .} in an effective 
way: the number z of a rational a/b can be computed from a and b, and vice versa. We fix 
such a numbering. 

Definition 2.3.2. A subset A of a numbered set O is recursively enumerable (r.e) if 

there is a r.e set C N such that A = {o„ : n G E}. j 

Uniformity for r.e subsets of O is defined as uniformity for r.e subsets of N. 



2.4. Computability of reals. The following notion was already introduced by Turing in 
[T36] . 

Definition 2.4.1. Let real number. We say that: 

• X is lower semi- computable if the set {g G (Q : g < a;} is r.e., 

• X is upper semi- computable if the set {g G (Q : g > x} is r.e., 

• X is computable if it is lower and upper semi-computable. j 

The following classical characterization may be more intuitive: a real number is com- 
putable if and only if there exists a recursive function (p computing a sequence of rational 
numbers converging exponentially fast to x, that is |g^(i) — x| < 2~*, for all i. We remark 
that as there exists subsets of integers which are recursively enumerable but not recursive 
Rog87| ), there also exists semi-computable numbers which are not computable. 



see 



In the following section we will see how these notions can be generalized to separable 
metric spaces, which inherit the computable structure of R via the metric. 
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2.5. Computable metric spaces. In this section we introduce the basic tools of com- 
putable analysis on metric spaces. Most of the results of this section and several of the 
following one have been already obtained by Weihrauch, Brattka, Presser and others in the 
framework of "Type-2 theory of Effectivity" , which is based in the notion of "representation" 
(infinite binary codes) of mathematical objects. A standard reference book on this approach 
to Computable Analysis is [WOO] . and a specific paper on computability of subsets of metric 
spaces is |BP03j . Our approach to Computable Analysis only uses the notion of recursive 
function (see subsection 12. 2p . It is intended to emphasize the fact that computability notions 
are just the "effective" versions of classical ones. In this way we obtain a theory syntactically 
familiar to most mathematicians and computability results can be proved in a transparent 
and compact way. 

A computable metric space is a metric space with a dense numbered set such that the 
distance on this set is algorithmically compatible with the numbering (distances between 
numbered points can be uniformly computed up to arbitrary precision). From this point of 
view the real line (with euclidean distance) has a natural structure of computable metric 
space, whit the rationals as a numbered set. 

Definition 2.5.1. A computable metric space (CMS) is a triple X = {X,d,S), where 

• (X, d) is a separable complete metric space, 

• S = (sj)jg]N is a dense subset of X (the numbered set of ideal points), 

• The real numbers {d{si, Sj))ij are all computable, uniformly in j 

Symbolic spaces, euclidean spaces, functions spaces and manifolds with a suitable metrics 
can be endowed with the structure of computable metric spaces. See for example |G93[ 
lHR09llGHR09b] . 

If {X, d, S) and {X', d', S') are two computable metric spaces, then the product {X x 
X',dx,S xS') with dx{{x,x'), {y,y')) = max((i(x, y) , d'{x', y') ) is a computable metric space. 

The numbered set of ideal points {si)i induces the numbered set of ideal balls B := 
{B{si,qj) : Si G S,qj G Q>o}- We denote by -B(i,j) the ideal ball B{si,qj). 

Let (X, d, S) be a computable metric space. The computable structure of X assures that 
the whole space can be "reached" using algorithmic means. Since ideal points (the finite 
objects of S) are dense, they can approximate any x at any finite precision. Then, every 
point X has a neighborhood basis consisting of ideal balls, denoted B{x) = {B G B : x E B} 
and called its ideal neighborhood basis. 

Definition 2.5.2 (Computable points). A point x E X is said to be computable if its ideal 
neighborhood basis B{x) is r.e. j 

Remark 2.5.1. As in the case of reals we have the following characterization: x is computable 
if and only if there is a (total) recursive function ip such that (i(s<^(i), x) < 2~\ _i 

Ideal balls are also useful to describe open sets. 

Definition 2.5.3 (Recursively open sets). We say that the set U C X is recursively open 

if there is some r.e set A of ideal balls such that U = [J^^^B. That is, if there is some r.e 
set E C N such that U = [J.^^ Bi. j 

We remark that the collection of r.e. open sets can be algorithmically enumerated. 
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Definition 2.5.4. Let {Un)n be a sequence of r.e. open sets. We say that the sequence is 
uniformly r.e. or that Un is r.e. open uniformly in n if there exists an r.e. set E' C 

such that for all n we have Un = IJigB where En = {i '■ {n, i) G E}. j 

Examples 2.5.1. 

(1) Let {Un)n be a sequence of open sets such that f/„ is uniformly recursively open. 
Then the union [J^ Un is a recursively open set. 

(2) The universal recursive function ip^ induces an enumeration of the collection U of 
all the recursively open sets. Indeed, define E := {(e, V9„((e, ra))) : e, n G N}. Then 
U = {U,:eeN} where U^ = U^eij, B^. 

(3) The numbered set U is closed under finite unions and finite intersections. Further- 
more, these operations are effective in the following sense: there exists recursive 
functions ip^ and ip'^ such that for all e, e' G N, ?7e U Ue' = f^ipLj(^e,e')) and the same 
holds for Equivalently: UeUUe' is recursively open uniformly in (e, e') (see |HR09] 
e.g.). 

Definition 2.5.5 (Computable functions). A function T : X ^ y is said to be computable 

if T^^{U^) is recursively open uniformly in e. j 

It follows that computable functions are continuous. Since we will work with functions 
which are not necessarily continuous everywhere, we shall consider functions which are com- 
putable on some subset of X. More precisely: 

Definition 2.5.6. A function T is said to be computable on C (C C X) if there is U^ 
recursively open uniformly in n such that 

T'\B:)nc = u^nc. 

The set C is called the domain of computability of T. j 

As an example we show that a monotone real function whose values over the rationals are 
computable, is computable everywhere. This Lemma will also be used later. 

Lemma 2.5.1. // / : [0,1] [0,1] is increasing and f{r) can be computed uniformly, for 
each rational r then f is computable. 

Proof. Let a, g G Q. We remark that f~^{{p, q)) = U/{a)>p,/(6)<g(a, b) this allows to find a r.e. 
cover of the interval f~^{{p,q))- The case of a general r.e. open set is straightforward. □ 

Definition 2.5.7 (Lower semi-computable functions). A function / : X IR is said to be 
lower semi- computable if f'^iqn, oo) is recursively open uniformly in n. _i 

It is known that there exists a recursive enumeration of all lower semi- computable functions 
{fi)i > 0. From the definition follows that lower semi-computable functions are lower semi- 
continuous. Lower semi- computability on D is defined as for computable functions. A 
function / is upper semi- computable if — / is lower semi-computable. It is easy to see 
that a real function / is computable if and only if it is upper and lower semi-computable. 

Given a probability measure /i, we say that a function is (lower semi-) computable 
almost everywhere if its domain of computability has /i-measure one. 

8 



2.6. Recursively compact sets: approximation from above. We will give some general 
results about solutions of equations concerning functions computable on some subset. As 
in many other mathematical situations, to prove the existence of certain solutions we are 
helped by a suitable notion of compactness. In order to the solution be computable, we will 
need a recursive version of compacity. Roughly, a compact set is recursively compact if the 
fact that it is covered by a finite collection of ideal balls can be tested algorithmically (for 
equivalence with the e-net approach see definition 12.7.11 and proposition 12.7.11 ) . This kind of 
notion and the related basic results are already present in the literature in various forms, or 
particular cases, we give a very compact self contained introduction based on the previously 
introduced notions. 

Definition 2.6.1. A set C X is recursively compact if it is compact and there is a 
recursive function : N — >• N such that ^{{ii, ■ ■ ■ , ip)) halts if and only if (Sj^, . . . , Bi^) is a 
covering of K. j 

Remark 2.6.1. Let Ut be the collection of r.e open sets (with its uniform enumeration). It 
is easy to see that a set K is recursively compact iS K O Ui is semi-decidable, uniformly in 
i. J 

Here are some basic properties of recursively compact sets: 

Proposition 2.6.1. Let K be a recursively compact subset of X . 

(1) A singleton {x} is recursively compact if and only if x is a computable point. 

(2) If K' is rec. compact then so is K U K' . 

(3) if U is recursively open, then K' = K \ U is rec compact. 

(4) The diameter of K is upper semi- computable. 

(5) The distance to K : dx{x) := mi{d{x,y) : y G K} is lower- computable 

(6) // / : X — > M lower- computable then so is inf^^ / 

(7) if f : X —* is upper- computable then so is sup^^ / 

Proof. (1) A point x is computable iff x G f/j is semi-decidable uniformly in i. {2) KUK' G U 
iS K C U and K' C U. (3) Remark that K\U CV K C U U V and U U V is 

recursively open uniformly in U and V. (4) diami^ = inf{g : 3s, K C B{s,q)}. (5) For x E X 
and g G Q define Ug^^ '■= {y '■ d{x,y) > q}, which is a constructive (in x) open set. Then 
dxi^) = sup{g : K C Ug^} is lower-computable. (6) inf^/ = sup{g : K C f~^[q,-\-oo)}. 
(7) snp^f = mf{q:KCf~\-oo,q)}. 

□ 

Remarks 2.6.1. 

(1) The arguments are uniform. In point 1) for instance, this means that there is an 
algorithm which takes a program computing x and outputs a program testifying the 
rec. compacity of {x}, and vice-versa. 

(2) When X itself is rec. compact, a subset K is rec. compact iff dx is lower-computable. 
Indeed, K = X\{x: dxi^x) > 0}. 

Corollary 2.6.1. // (-ft'j)jgN are uniformly recursively compact sets, then so is f]-^^Ki. 

Proof. The complements of recursively compact sets are r.e open. Then by proposition 12. 6. 1| 

part ([2]) the set HieN -^i = -^o \ (Ui>o ^i) recursively compact. □ 
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It is important to remark that a recursively compact set needs not contain computable 
points. This will be used in section [51 



Proposition 2.6.2. There exists a nonempty recursively compact set K C [0, 1] containing 
no computable points. 

Proof. Let be an enumeration of all the rational intervals and e > be a rational number. 
Put E = {i >1 : (Pi{i) halts and < e2~*}. E is a. i.e. subset of N. Let U = [J.^^Ii: 

X{U) < ^jg£; e2~* < e. Let x G [0, 1] be a computable real number. There is a total recursive 
function ip^ such that < e2~" and x G I<p^{n) for all n. Hence i E E, so x eU. Hence 

U contains all computable points. As [0, 1] is recursively compact, so is if = [0, 1] \ t/. □ 

Now we start to show that many statements about topology and calculus on metric spaces 
can be easily translated to the computable setting: the first one says that the image of a 
recursively compact is still a recursively compact. 

Proposition 2.6.3 (Stability by computable functions). Let f : K (1 X —>■ Y be a com- 
putable function defined on a recursively compact set K . Then f{K) is recursively compact. 



Remark that the argument is uniform: if {Ki)i(z-^ is a sequence of uniformly recursively 
compact subsets of X on which / is defined, then {f{Ki))i(z^ is a sequence of uniformly 
recursively compact subsets of Y . We will say that f{K) is recursively compact uniformly 
in K. 

As a first simple example of application, we observe that in some cases the global attractor 
of a (computable) dynamical system can be approximated by an algorithm to any given 
accuracy. 

Corollary 2.6.2. Let X be a recursively compact computable metric space and T a com- 
putable dynamics on it. Then the set: 



We remark that these and other frameworks of "exact computability and rigorous approx- 
imation" have been previously used to study the computability of several similar objects 
such as Juha or Mandelbrot sets f [H05l IBY061 IBBY061 IBBYOTj . |Del97 ]). or the existence 
and some basic properties of Lorentz attractor ( |Tuc99] ) . 

Here is a computable version of Heine's theorem. 

Definition 2.6.2. A function f : X ^Y between metric spaces is recursively uniformly 
continuous if there is a recursive 5 : Q — ^ (Q such that for all e > 0, (5(e) > and Vx G X, 



Proof. Indeed, f{K)CU ^ Is 
total recursive function, f{K) C U, 




A := Pi T"(X) 




Proof. By proposition 12.6.31 and corollary 12.6.11 



□ 




(2.1) 



J 
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Proposition 2.6.4. Let X and Y be two computable metric spaces. Let K C X be recur- 
sively compact and f : K ^ Y be a computable function. Then f is recursively uniformly 
continuous. 

Proof. First, x i^' is a recursively compact subset of X x X. For eacli rational number 
e > 0, define f/(e) = {{x,x') e : d{f{x)J{x')) < e] and K{t) = KxK\ U{e): they are 
respectively recursively open and recursively compact, uniformly in e. Hence, the function 
S{e) := mi{d{x,y) : {x,y) G K{e)} is lower semi-computable (proposition [6]) . 

Now, / is uniformly continuous if and only if (5(e) > for each e > 0. By the classical 
Heine's theorem, this is the case, so by lower semi-computability of S{e), one can compute 
from e some positive 6 < S{e). □ 

Theorem 2.6.1. Let K be a recursively compact subset of X and f : K ^ H be a computable 
function. Then every isolated zero of f is computable. 

Proof. Let Xq be an isolated zero of /. Let s,r be an ideal point and a positive rational 
number such that x G B{s,r) and the only zero of / lying in B{s,r) is Xq. The set N = 
{x : f{x) ^ 0} U {x : d{x, s) > r} is recursively open in K (that is, N H K = U H K with 
U recursively open), so {xq} = K\ N = K\ U is recursively compact by proposition 12.6. II 
Hence, Xq is a computable point. □ 

Remark 2.6.2. Observe that the argument is uniform in / and an ideal ball isolating the 
zero. In particular, there is an algorithm which takes a finite description of / and the ball 
and outputs his zero if it is unique. j 

Corollary 2.6.3. Let K be a recursively compact subset of X and f : K X be a com- 
putable function. Then every isolated fixed point of f is computable. 

Proof. Apply the preceding theorem to the function (7 : X — > R defined by g{x) = d{x, f{x)). 

□ 

2.7. Recursively precompact. In this subsection we prove the equivalence between the 
notion of recursive compactness given above and another natural approach (which will be 
used later) to recursive compactness, where it is supposed the existence of an algorithm to 
construct e-nets. 

Definition 2.7.1. A CMS is recursively precompact if there is a total recursive function 
99 : N — >• N such that for all n, ip{n) computes a 2~"-net: that is (p{n) = {ii, . . . ,ip) where 
(si,, . . . , Si^) is a 2""-net. j 

Here is a computable version of a classical theorem: 

Proposition 2.7.1. Let X be a CMS. X is recursively compact if and only if it is complete 
and recursively precompact. 

Proof. If X is recursively compact then we define the following algorithm: it takes n as input, 
then enumerates all the {ii, . . . ,ip), and tests if (-B(si^, 2""), . . . , B^Si^, 2~")) is a covering 
of X (this is possible by recursive compacity). As X is compact, hence precompact, such 
a covering exists and will be eventually enumerated: output it. The algorithm makes X 
recursively precompact. 

Suppose that X is complete and recursively precompact. Let {B{si, qi), . . . , B{sk, qu)) be 
ideal balls: we claim that {B{si, gi), . . . , B{sk, qk)) covers X if and only if there exists n such 
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that each point s of the 2~"-net given by recursive precompactness hes in a ball B{si, qi) 
satisfying (i(s, Si) + 2~" < g^. The procedure which enumerates all the n and semi-decides 
this halts if and only if the initial sequence of balls covers X. We leave the proof of the 
claim to the reader (take n such that 2~" is less than the Lebesgue number of the finite 



The following observation is also worth noticing. 

Proposition 2.7.2. Let X he a computable metric space. If X (as a subset of X) is recur- 
sively compact, then the set C{X) of continuous functions from X to 'K with the distance 
induced by the uniform norm is a computable metric space. 

The function eval : C{X) x X — > R mapping (/, x) to f{x) is computable. 

Let Y be a computable metric space: for every computable function f : Y x X H, the 
function Y — > C{X) mapping y to fy : x \—>- f{y, x) is computable. 

2.8. Recursively closed sets: approximable from below. From the computability 
viewpoint, the properties of recursively closed sets are, in a sense, complementary to those 
of recursively compact sets. 

Definition 2.8.1. A closed set F is recursively closed if the set {B{s, r) : B{s, r)nF ^ 0} 
is r.e. J 

A closed set F is recursively closed if F fl f/ is semi-decidable for r.e open sets U . It is easy 
to see that the union of two recursively closed sets is also recursively closed. The closure of 
any recursively open set is recursively closed: -B fl f/ 7^ <^=^ 3s G -B fl U . 

The following proposition will be used later. 

Proposition 2.8.1. Let F be a recursively closed subset of X. Then there exists a sequence 
of uniformly computable points Xi & F which is dense in F . 

Proof. Since {ra G N : -B„ = g„) fl F 7^ 0} is r.e, given some ideal ball B = B{s, q) 

intersecting F, the set {n G N : i?„ C 5, g„ < 2~", S„ fl F 7^ 0} is also r.e. Then we 
can effectively construct an exponentially decreasing sequence of ideal balls intersecting F. 
Hence {x} = HkBh is a computable point lying in F. 



We remark that by this. Proposition 12.6.21 shows a recursive compact which is not recur- 
sively closed. For the sake of completeness, let us state some useful simple properties. 

Proposition 2.8.2. Let F be a recursively closed subset of X . Then: 

(1) The diameter of F is lower semi- computable, uniformly in F. 

(2) If f : F H is lower semi- computable, then so is suppf. 

(3) If f : F ^ H is upper semi-computable, then so is inip f. 

Proof. (1) Let C{s,q) be the complement of the closed ball B{s,q), that is C{s,r) = {x : 
d{x,s) > q}: this is a recursively open set, uniformly in s,q. Then diamF = sup{g : 



3s, C{s, r) n F ^ 0}. (2) sup^ / = sup{g : f-\q, +00) H F ^ 0}. (3) Apply (2) to -/. □ 



Corollary 2.8.1. Let K be recursively closed and recursively compact subset of X . If f : 
K is a computable function, then so are mix f and supj^ /. 



covermg ) . 



□ 



□ 
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3. Computable measures 



Let us consider the space PM{X) of Borel probability measures over X. We recall that 
PM{X) can be seen as the dual of the space Co{X) of continuous functions with compact 
support over X and recall the notion of weak convergence of measures: 

Definition 3.0.2. /i„ is said to be weakly convergent to /i if J /d/i„ — > j f dfi for each 
/ G Co{X). 

Let us introduce the Wasserstein-Kantorovich distance between measures. Let /i^ and /ig 
be two probability measures on X and consider: 

Wi{fii,fi2) = sup fdfi^- fdfi 
/ei-Lip{x) J J 

where l-Lip(X) is the space of 1-Lipschitz functions on X. We remark that since adding 
a constant to the test function / does not change the above difference / / dfii — J f d/Xj then 
the supremum can be taken over the set of 1-Lipschitz functions mapping a distinguished 
ideal point Sq to 0. The distance Wi has moreover the following useful properties which will 
be used in the following 

Proposition 3.0.3 f jXGS] Prop 7.1.5). 

(1) Wi is a distance and if X is bounded, separable and complete, then PM{X) with this 
distance is a separable and complete metric space. 

(2) // X is bounded, a sequence is convergent for the Wi metrics if and only if it is 
convergent for the weak topology. 

(3) If X is compact PM{X) is compact with this topology. 

Item (1) has an effective version: PM{X) inherits the computable metric structure of X. 
Indeed, given the set Sx of ideal points of X we can naturally define a set of ideal points 
SpM(x) in PM{X) by considering finite rational convex combinations of the Dirac measures 
5x supported on ideal points x & Sx- This is a dense subset of PM{X). The proof of the 
following proposition can be found in ( |HR09] ) 

Proposition 3.0.4. If X bounded then PM{X) with the Wi distance (and SpM{x) as a set 
of ideal points) is a computable metric space. 

A measure /i is then computable if there is a fast sequence G SpM{x) converging to 
/i (see remark r2. 5. ip in the Wi metric (and hence for the weak convergence). 
Now, point (3) of proposition I3.U.3I also has an effective version: 

Lemma 3.0.1. If X is a recursively precompact metric space, then PM{X) with the Wi 
distance is a recursively precompact metric space. 

Proof. We will show how to effectively find an r— net for each r of the form r = ;i,n G N. 
Let us consider the set S*,- = {^, < A; < n} subdividing the unit intervals in equal segments. 
Let us also consider an r-net Nr = {xi, ...Xm} constructed by recursive compactness of X. 
Now let us consider the set of measures with support in A^^ given by 

Tj- = {ki5xi ~\~ ... + kfnSxm s.t. ki G Sj- , ki + ... + km — 1}- 

This is a 2r net in PM{X). To see this let us consider a measure /i on X and a ball B{xi, r) 
centered in xi G X. Let us consider the measure fii defined by 
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for each measurable set A G X. The measure is obtained transporting the mass contained 
in the ball B{xi,r) to its center. Then Wi{fii,fi) < rfi{B{xi,r)). Let us now consider the 
sequence of measures fii, where fii is as before and the other ones are given by 

at the end /i„ is a measure with support in Nr and by the triangle inequality Wi{fij^, /i) < r. 

Now /i^ has the same support as the measures in and there is G such that 
I / f dfi^ — / /dz/| < r for each / G l-Lip(X), hence Wi{fj.j^, u) < r and then Wi{fi, u) < 2r 
and this proves the statement. □ 

We now use the recursive enumeration of lower semi- computable functions > to 
characterize computability on PM{X) (see |HR09j corollary 4.3.1): 

Lemma 3.0.2. Let X be a bounded computable metric space and S be any subset of PM{X), 
then: 

(1) /i G PM{X) is computable iff the function fi J fid^ is lower semi-computable, 
uniformly in i, 

(2) L : PM{X) — > PM{X) is computable on S iff the function ^ ^ J fi dL{fi) is lower 
semi-computable on S, uniformly in i. 

This gives: 

Lemma 3.0.3. If Qi : X —>■ R"*" is a uniform sequence of functions which are lower semi- 
computable on X\D, then ^ Jgidfi is lower semi- computable on 

PMniX) := {fi : fi{D) = 0} (3.1) 

uniformly in i. 

Proof. For each i, one can construct a lower semi-computable function Qi satisfying gi = on 
X\D (see |HR09] . subsection 3.1). Since the function yU t— > jgidfi is lower semi-computable, 
uniformly in i and fi{D) = 0, we have that on PM£){X) it coincides with fi Jgi d/i, which 
is then lower semi-computable on PMij{X), uniformly in i. □ 

An interesting remark about computable measures is that they must have computable 
points in the support. This will be used in section [STTl 

Proposition 3.0.5. If ^ is a computable probability measure, then there exists computable 
points in the support of fi. 

Proof. The sequence of functions fi := 1^^ (the indicator functions of ideal balls) are uni- 
formly lower semi-computable. By lemma 13.0.21 the numbers J fidfi = fi{Bi) are uniformly 
lower semi-computable. Hence, the set {Bi : n{Bi) > 0} is recursively enumerable. In other 
words, the support of /i is a recursively closed set. Proposition 12.8.11 allows to conclude. □ 
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4. Dynamical systems, statistical behavior, invariant measures 

Let X be a metric space, let T : X i-^ X be a Borel measurable map. Let fi be an invariant 
measure. A set a A is called T-invariant if T~^{A) = A{modO). The system {T,fi) is said 
to be ergodic if each T-invariant set has total or null measure. In such systems the famous 
Birkhoff ergodic theorem says that time averages computed along fi typical orbits coincides 
with space average with respect to fi. More precisely, for any / G L^{X) it holds 

liin = [ fdfi, 

n— >oo n J 

for /i almost each x, where Sl = f + foT + ... + foT' 

This shows that in an ergodic system, the statistical behavior of observables, under typical 
realizations of the system is given by the average of the observable made with the invariant 
measure. 

In case X is a manifold (possibly with boundary). We say that a point x belong to the 
basin of an invariant measure fi if Equation 14.11 holds at x for each continuous / (the average 
on the X orbit represent the average under the measure). An SRB measure is an invariant 
measure having a positive Lebesgue measure basin (for more details and a general survey 
see |Y02] l 

In the applied literature the most common method to simulate or understand the above 
statistical behaviors is to compute and study some trajectory. This method has three main 
theoretical problems which motivates the search of another approach: 

• numerical error, 

• tipicality of the sample, 

• how many sample points are necessary? 

the first (and widely known) problem is the amplification of the numerical error (if the 
system is sensitive to initial conditions as most interesting systems are). Here the shadowing 
results are often invoked to justify the correctness of simulations, but rigorous results are 
proved only for a small class of systems (see e.g. |PalOO] ) and moreover the mere existence of 
a shadowing orbit does not say anything about its typicality (see e.g. |Bla89t IBla94j for a 
further discussion on numerical errors). 

The second problem is indeed that this method should compute, in order to be useful, a 
trajectory which shows the "typical" behavior of the system: a behavior which take place 
with large or full probability in some sense. The main problem here is the fact that the 
set of initial conditions the computer has access to, being countable, has probability zero. 
Hence, there is no guarantee that what we see on the screen is typical in some sense. On 
the contrary, in a chaotic system, typical orbits are far from being describable by a finite 
program. It is true for example that in an ergodic system having positive entropy h a typical 
n step orbit segment needs approximatively a program which is hn bits long to be described 
(up some approximation e, see e.g. |B83] for the original result or |GaOO] and [GHROQc] for 
a version in the framework of computable analysis). We remark, however, that if one looks 
for points which behave as typical for Birkhoff averages (hence they behave as typical for 
some given particular aspect) there are some rigorous results partly supporting this way to 
proceed: in several classes of systems there are computable initial conditions which behave 
as typical with respect to Birkhoff averages (see |GHR09b] for a precise result). 
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The third problem however remains. Even if you find a program describing a typical orbit 
of the system: how many iterations should be considered to be near to the limit behavior, 
so that the orbit represents the invariant measure up to a certain approximation? although 
this problem can be approached rigorously in some cases (see |CCS] e.g.) we will not adopt 
this point of view. We will study the system's statistical behavior by directly computing the 
invariant measure as fixed points of a certain transfer operator. 

4.1. The transfer operator. A function T between metric spaces naturally induces a 
function between probability measure spaces. This function Lt is linear and is called transfer 
operator (associated to T) . Measures which are invariant for T Invariant measures are fixed 
points of Lt- 

Let us consider a computable metric space X endowed with a Borel probability measure 
/i and with a dynamics defined by a measure-preserving function T : X ^ X. Let us also 
consider the space PM{X) of Borel probability measures on X. 

Let us define the function Lt '■ PM{X) PM{X) by duality in the following way: if 
fi e PM{X) then LT(/i) is such that 

y"/dLT(/i) = jfoTdfx 

for each / G Cq{X). In next sections, invariant measures will be found as solutions of the 
equation Wi{fi, L{fi)) = 0. To apply Theorem 12.6.11 and Corollary 12.6.31 to this equation we 
need that L is computable. We remark that if T is not continuous then L is not necessarily 
continuous (this can be realized by applying L to some delta measure placed near a discon- 
tinuity point) hence not computable. Still, we have that L is continuous (and its modulus of 
continuity is computable) at all measures /i which are "far enough" from the discontinuity 
set D. This is technically expressed by the condition /i(-D) = 0. 

We remark that with the general tools introduced before, the proof is immediate. 

Theorem 4.1.1. Let X be a computable metric space and T : X ^ X be a function which 
is computable on X \D. Then Lx is computable on the set of measures 

PMd{X) := {/i G PM{X) : /i(D) = 0}. (4.2) 

Proof. Note that if / is lower semi-computable, then foT is lower semi-computable on X\D. 
The result then follows from lemmas 13.0.31 and I3.U.2I □ 

In particular, if T is computable on the whole space X then L is computable on all 
PM{X). 

4.2. Computing invariant "regular" measures. The above tools allow to ensure the 
computability of Lt on a large class of measures. This will allow to apply Corollary 12.6.31 
and see an invariant measure as a fixed point. 

Theorem 4.2.1. Let X be a computable metric space and T be a function which is com- 
putable on X \D. Suppose there is a recursively compact set of probability measures V C 
M{X) such that for every fi & V , yu(-D) = holds. Then every invariant measure isolated in 
V is computable. 
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Proof. By Theorem 14.1.11 Lt is computable on V. Since V is recursively compact, theorem 
12.6.11 allows to compute any invariant isolated measure in as a solution of the equation 
Lxifi) = /i. □ 

Remark 4.2.1. This theorem is uniform: there is an algorithm which takes as inputs finite 
descriptions of T, V and an ideal ball in M{X) which isolates an invariant measure /x, and 
outputs a finite description of fi (see the above proof and Remark [2.6.21) . j 

A trivial consequence of Theorem 14.2.11 is the following: 

Corollary 4.2.1. // a computable system as above is uniquely ergodic and its invariant 
measure fi satisfy ii{D) = 0, then it is a computable measure. 

The main problem in the application of theorem 14.2.11 is the requirement that the invariant 
measure we are trying to compute, should be isolated inV. In general the space of invariant 
measures in a given dynamical system could be very large (an infinite dimensional convex in 
PM{X) ) to isolate a particular measure we can restrict and consider a subclass of "regular" 
measures. 

Let us consider the following seminorm: 

II II KB{x,r)) 

\\^'L= sup — - — . 

xex,r>o r" 

Proposition 4.2.1. If X is recursively compact then 

Va,K = {/i e PM{X) ■ < K} (4.3) 

is recursively compact. 

Proof. U = {fi E PM{X) : > K} is recursively open. Indeed, > K iff there 

exists s, r G iS X Q for which fi{B{s,r)) > qr". As fi ^ fi{B{s,r)) is lower semi-computable 
uniformly in s,r, the sets Us^r '■= {f^ '■ fJ'{B{s,r)) > Kr"'} are uniformly recursively open 
subsets of PM{X). Hence, U = Us^rUs,r is recursively open. 

Now, Va^K = PM{X) \ U . As PM{X) is recursively compact (see Lemma [3.0.11) and U is 
recursively open, then proposition 12.6.1] part (3) allows to conclude. □ 

In theorem 14.2.11 we require that /i(-D) = holds. This is automatically true in many 
examples when the measure is regular and the set D is reasonably small. 

Proposition 4.2.2. Let X be recursively compact and T be computable on X \ D , with 
dimuiD) < oo. Then any invariant measure isolated in Va^K with a > dimniD) is com- 
putable. 

Proof. Let us first prove that IJ^{D) = for all /i G Va^K- For all e > 0, there is a covering 
{B{xi,ri))i of D satisfying ^.rf < e. Hence fi{D) < '^^ /iS(xi, n) < 2"i^ ^.rf < 2-ife. 
As this is true for each e > 0, /i(-D) = 0. 

The result then follows from the fact that Va,K is recursively compact and Theorem 14.2.11 

□ 

Remark 4.2.2. Once again, this is uniform in T, a, K. j 

The above general proposition allows to obtain as a corollary the computability of many 
absolutely continuous invariant measures. For the sake of simplicity, let us consider maps 
on the interval. 
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Proposition 4.2.3. If X = [0, 1], T is computable on X\D, with dimH{D) < 1 and {X,T) 
has an unique a.c.i.m. fi with bounded density, then ^ is computable. 

Proof. The result follows from the above proposition 14.2.21 and the fact that if is absolutely 
continuous and the density of is / e Li[0, 1] then = esssup{f). We have to check 
that there could be not other measures having a finite 1 norm and not being absolutely 
continuous. 

If we suppose that = / is finite, then /i is absolutely continuous, with bounded density 
f < I. Indeed, let us consider the conditional expectation £'[/i|J"] of /i to the dyadic n-th 
grid J" = {[fc2-", {k + 1)2""), < A; < 2"}. 

If = / a fortiori implies < -E'[/i|/"] < / a.e.. By the first Doobs martingale conver- 
gence it follows that has an a.e. pointwise limit / and / < / a.e.. Since / is bounded 
then it is a density for /x. 

□ 

(i-dimensional submanifolds of R" can naturally be endowed with a natural structure of 
computable metric spaces ( see |GHR09b] ). Considering a dyadic grid on and chart 
diffeomorphisms it is straightforward to prove, in the same way as before 

Corollary 4.2.2. Let X be a recursively compact d dimensional submanifold o/R" (with 
or without boundary). If T is computable on X\D, with dimn^D) < d and {X,T) has a 
unique a.c.i.m. fi with bounded density, then /i is computable. 

As it is well known, interesting examples of systems having an unique a.c.i.m. (with 
bounded density as required) are topologically transitive piecewise expanding maps on the 
interval or expanding maps on manifolds (see |V97] for precise definitions). Provided that 
the dynamics is computable we then have by the above propositions that the a.c.i.m. is 
computable too. 

4.3. Unbounded densities. The above results ensure computability of measures having 
an a.c.i.m. with bounded density. If we are interested in situations where the density is 
unbounded, we can consider a new norm, "killing" singularities. 
Let us hence consider a computable function / : X — > R and 

_ f{x)fx{B{x,r)) 

^' x€X,r>0 

Propositions 14.2.11 and 14.2.21 also hold for this norm. If / is such that /(x) = when 
lim^^o ^^^IZ'^^'' = oo this can let the norm to be finite when the density diverges. 

As an example, where this can be applied, let us consider the Manneville Pomeau maps 
on the unit interval. These are maps of the type x — x + x^(mod 1). When 1 < z < 2 
the dynamics has an unique a.c.i.m. /i^ having density ez{x) which diverges in the origin as 
ez{x) X x~^~^^ and it is bounded elsewhere (see |I03] section 10 and |V97] section 3 e.g.). If 
we consider the norm \\-\\fi with f{x) = x"^ we have that H/U^Uji is finite for each such z. By 
this it follows that the measure /i^ is computable. 

5. Computable systems having not computable invariant measures 



We have seen that the technique presented above proves the computability of many a.c.i.m. 
which are also SRB measures. As we have seen in the introduction, with other techniques 
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it is possible to prove the comput ability of other SRB measures (axiom A systems e.g., 
see |GHR09b] ). This raises naturally the following question: a computable systems does 
necessarily have a computable invariant measure? what about ergodic SRB measures? 

The following is an easy example showing that this is not true in general even in quite 
regular systems, hence the whole question of computing invariant measures has some subtlety. 

Let us consider a system on the unit interval given as follows. Let r G (0, 1) be a lower 
semi-computable real number which is not computable. There is a computable sequence 
of rational numbers such that supj Tj = r. For each i, define Ti{x) = max(a;, r,) and 
T(x) = Xlj>i2~*2^i- The functions Tj are uniformly computable so T is also computable. 




Figure 1 . The map T. 
Now, T is non-decreasing, and T{x) > x if and only if x < r. 

The system ([0, 1],T) is hence a computable dynamical system. This system has a SRB 
ergodic invariant measure which is 6r, the Dirac measure placed on r. The measure is SRB 
because r attracts all the interval at its left. Since r is not computable then 6r is not 
computable. We remark that coherently with the previous theorems 6r is not isolated. 

We remark that by a simple dichotomy argument we can prove that a computable function 
from [0, 1] to itself must have a computable fixed point. Hence it is not possible to construct 
a system over the interval having no computable invariant measure (we always have the 
6 over the fixed point). With some more work we will see that such an example can be 
constructed on the circle. 

5.1. A computable system having no computable invariant measure. We go further 
and exhibit a computable dynamical system on a compact space which has no computable 
invariant probability measure. 

We consider the unit circle S, identified with M/Z. It naturally has a computable metric 
structure inherited from that of M. 

On S, there is a computable map with no computable invariant probability measure. We 
construct such a map T : [0, 1] M satisfying T(l) = T(0) + 1, and consider its quotient on 
the unit circle. 

From proposition 12.6.21 we know that there is a non-empty recursively compact set K 
containing no computable point. Let U = (0, 1) \ K: this is a r.e. open set, so there are 
computable sequences aj,6j {i > 1) such that < Oj < 6,, < 1 and U = |Jj('^«5^i)- Let us 
define non-decreasing, uniformly computable functions /, : [0, 1] — > [0, 1] such that fi{x) > x 
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if a; e {ai,bi) and fi{x) 
fi{x)^h on [^Ml 



Figure 2. The map /j. 

As neither nor 1 belongs to K, there is a rational number e > such that K C. [e, 1 — e] . 
Let us define / : [0, 1] — > IR by f{x) — x on [e, 1 — e], f{x) — 2x — {1 — e) on [1 — e, 1] and 

f{x) = e on [0,e]. 

We then define the map T : [0, 1] — ^ M by T{x) = ^ + Xli>2 '^^^fi- ^ computable and 
non-decreasing, and T{x) > a; if and only ii x E [0,1]\K. As T(l) = 1 + | = 1 + T(0), we 
can take the quotient of T modulo 1. 




Figure 3. The map T. 

Proposition 5.1.1. W — U U [0, e) U {1 — e, 1] is a strictly invariant set: T'^W — W. 

Proof. lix^W then T(x) ^x^W. 

U x E W then T{x) G W. Indeed, if T{x) ^ W, T{x) is a fixed point so T is constant on 
[x, T(a;)] (T is non-decreasing). Let q be any rational number in {x,T{x)): T[x) = T[q) is 
then computable, but does not belong to W: impossible. □ 

Proposition 5.1.2. The map T is computable hut has no computable invariant probability 
measure. 

Let X e [0,1]: the trajectory of x is "non-decreasing" and converges to the first point 
above x which is not in U , inf([a;, 1] \ U) or to min(i^r) if x > sup(i^'). More precisely, there 
are two cases: (i) \i x then x is a fixed point (unstable on the right), (ii) if x e C/ then 
the trajectory of x converges to a lower semi-computable fixed point (non-computable, as it 
does not belong to U). 

Lemma 5.1.1. Let /i be an invariant probability measure: then n{K'^) — 0. 

Proof. Obviously ;u(0) = because is not periodic. Let (a, h) — (oj, h,i) be an interval from 
the description of U . Since T"-{a) and T"(6) tends to some non computable a (and then are 
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= X otherwise. For instance, fi{x) — 2x — ai on [a^, 2i±2i] and 




not stationary, as they are computable), the interval (a, b) is wandering. Hence, by Poincare 
recurrence theorem it has null measure. □ 



Proof, (of proposition 15.1.21) We can conclude: let /i be a computable invariant probability 
measure: by the above lemma its support is then included in the complement of W. But 
the support of a computable probability measure always contains computable points (see 
proposition 13. 0.Sp : contradiction. □ 

Actually, the set of invariant measures is exactly the set of measures which give null 
weight to W. It is easy to see that in the above system the set of invariant measures is a 
convex recursive compact set. Indeed, the function fi — > fi{W) is lower semi-computable, so 
{ij, : fJ,{W) > 0} is a recursive open set. Its complement is then a recursive compact set, 
as the whole space of probability measures is a recursive compact set. The above example 
hence shows an example of a convex, and recursive compact set whose extremal points are 
not computable. 

We end remarking that with a different construction of the various fi it is possible to give 
also a smooth system having the same properties as the examples in this section. 
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